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Abstract 



We consider the zeroes of the random Gaussian entire function 



z 

fc=0 A 



(Co ? 5 • • • are Gaussian i.i.d. complex random variables) and show that their basins 
under the gradient flow of the random potential U(z) = log \f(z)\ — \\z\ 2 partition 
the complex plane into domains of equal area. 

We find three characteristic exponents 1, |, and 4 of this random partition: the 
probability that the diameter of a particular basin is greater than R is exponentially 
small in R; the probability that a given point z lies at a distance larger than R from 
the zero it is attracted to decays as e~ R8/ °; and the probability that, after throwing 
away 1% of the area of the basin, its diameter is still larger than R decays as e~ R . 

We also introduce a combinatorial procedure that modifies a small portion of each 
basin in such a way that the probability that the diameter of a particular modified 
basin is greater than R decays as e ~ cR4 ( l °sR) 3/2 _ 
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1 Introduction and main results 



Let Z be a random point process in IR d with the distribution invariant with respect 
to the isometries of M. d . Suppose that Z has intensity 1; that is, the mean number 
of points of Z per unit volume equals 1. The transportation (a.k.a. "matching", 
"allocation", "marriage", etc.) of the Lebesgue measure to Z is a (random) 
measurable map T : M, d — > Z that pushes forward the Lebesgue measure to the 
counting measure nz = S a of the set Z ( S a is the unit mass at a). In other words, 

the whole space M. d is split into disjoint random sets B(a) of the Lebesgue measure 1 
indexed by a e -Z. Because of the invariance of the process Z, it is natural to assume 
that the transportation T has an invariant distribution; i.e., that the distribution 
of the vector T(x) — x does not depend on x. The better T is localized, the more 
uniformly the process Z is spread over M. d . Thus it is interesting to know the optimal 
rate of decay of the probability tails F {\T(x) — x\ > R} as i? ^ oo. A constructive 
counterpart is to find an explicit and well-localized way to transport the Lebesgue 
measure to the point process Z. 

The transportation of the Lebesgue measure to the Poisson process in M. d was 
recently developed by Hoffman, Holroyd and Peres [HI E] (a finite volume version 
was studied earlier by Ajtai, Komlos and Tusnady [T], Leighton and Shor [13], and 
Talagrand [17]) . In this paper, we consider the random zero point set Zf = / _1 (0) of 
a Gaussian entire function / in C and study the transportation of the two-dimensional 
Lebesgue measure m 2 to Zf. 

Let 



fc>0 v 

where are independent standard complex Gaussian random variables (i.e., the 
density of £/. on the complex plane C is )• We shall call such a random function 

a Gaussian Entire Function (G.E.F.). 

The (random) zero set Zf of this function is known as "flat chaotic analytic zero 
points" P, El H2] • It is distinguished by the invariance of its distribution with respect 
to the isometries of C; i.e., rotations and translations, see [151 Part I] for details and 
references. Note that the intensity of the zero process Zj equals -. In [T5l part II], the 
question about the existence of a well-localized transportation of the area measure 
to the zero set of the Gaussian Entire Function in C was studied. Using the Hall 
matching lemma and some potential theory, the authors of [T3] proved the existence 
of a transportation with sub-Gaussian decay of the tail probability. Unfortunately, 
the proof one obtains on this way is a pure existence proof giving no idea of what the 
transportation in question looks like. 
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The aim of this paper is to carry out another approach that was suggested but not 
followed in (TSJ part II] , namely, the transportation by the gradient flow of a random 
potential. The main advantage of this approach is that it provides a quite natural 
and explicit construction for the desired transportation. 

Let U(z) = log \f(z)\ — ^\z\ 2 be the random potential corresponding to the G.E.F. 
/. The distribution of U is also invariant with respect to the isometries of the complex 
plane, see [To"} P ar t I] or Section 2.2 below. We shall call any integral curve of the 
differential equation 

dZ , . 

n = -w 

a gradient curve of the random potential U. 

We orient the gradient curves in the direction of decay of U (this is the reason 
for our choice of the minus sign in the differential equation above). If z £ Zf, and 
VU(z) 7^ 0, by Y z we denote the (unique) gradient curve that passes through the 
point z. 

Definition 1.1 (the basin) Let a be a zero of the G.E.F. f. The basin of a is the 
set 

B(a) = {z G C: VU(z) ^ 0, and T z terminates at a} . 

The picture below may help the reader to visualize this definition. It shows the 
random zeroes and the trajectories of various points under the gradient flow. 




Figure 1: The basins B(a) and trajectories of the gradient field 

Clearly, each basin B(a) is a connected open set, and B(a!) H B{a") = if a' and 
a" are two different zeroes of /. If the basin B(a) is bounded and the boundary of 
B(a) is nice, then |^ = on dB(a) and therefore, applying the Green formula and 
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observing that the distributional Laplacian of U equals AU = 27r5^ oeZ <J„ — 2m2, 
one gets 

TT 2nJjB(a) 2nJdB(a)dn 

i.e., rri2B(a) = ir. 

Now we are ready to formulate our main results: 

Theorem 1.2 (partition) Almost surely, each basin is bounded by finitely many 
smooth gradient curves (and, thereby, has area n), and 

C = |J B ^ 

up to a set of measure (more precisely, up to countably many smooth boundary 
curves) . 

Consider the random set 

S = (J dB{a) ■ 

aeZ f 

that is, the union of all "singular" gradient curves; i.e., the curves that do not ter- 
minate at Zf. Due to the translation invariance of the random potential U, the 
probability P { z G S } does not depend on the choice of the point z G C, hence 
vanishes: 

P{0gS} = - // F{zeS} dm 2 (z) = [ m 9 (S n D) dF = 

(here f2 is the probability space and D is the unit disk). Thus, almost surely, any 
given point z G C belongs to some basin. 

By B z we denote the basin that contains the point z. By diam(A) we denote the 
diameter of a set A C C. We denote by C and c absolute (numerical) constants that 
may change from one line to another. 

Theorem 1.3 (diameter of the basin) For any point zgC and any R ^ 1, 

ce -CR(io gR) ^ ^ p{ diam ( jB2 ) > R y <; ce-c^iog*)^ 

The proof of Theorem 11.31 relies on the following auxiliary theorem. Let Q(w, s) 
be the square centered at w with side length 2s and let dQ(w, s) be its boundary. 
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Theorem 1.4 (long gradient curve) Let R ^ 1. T/ie probability of the event 
that there exists a gradient curve joining dQ(0, R) with dQ(0, 2R) does not exceed 

(j e -cB.(\og _R) 3 / 2 

The proof of this theorem is, unfortunately, quite involved. For a weaker upper 
bound Ce~ cR ^ logR that has a simpler proof, see the first version of this work posted 
in the arxiv [T3j. The approach in [H] may be more suitable for extensions to point 
processes of different nature: recently, using a similar approach, Chatterjee, Peled, 
Peres, and Romik found counterparts of Theorems 1 1 . 3 1 and 1 1 . 41 for the Poisson process 
in M. d with d ^ 3 [3]. It might be helpful for the reader to look at [H] prior to reading 
the proof of the long gradient curve theorem given here. 

Let a z be the random zero whose basin contains a given point z G C. In other 
words, the gradient curve T z terminates at a z . It appears that the probability P{ \z — 
a z \ > R } is much smaller than the probability P { dian^i^) > R } : 

Theorem 1.5 (distance to the sink) For any point zeC and any R ^ 1, 

ce- CR8/5 ^ ¥{\z-a z \>R} ^Ce~ cR8/5 . 

This is related to long, thin "tentacles" seen on the picture around some basins. 
They increase the typical diameter of basins though the probability that a given point 
z lies in such a tentacle is very small. 

Let D(w,r) be the disk of radius r centered at w. 

Theorem 1.6 (diameter of the core) For any z G C ; any e > 0, and any R ^ 1, 

c(e)e- c ^ R4 ^¥{m 2 (B z \ D(a z , R))^e} ^ C(e)e- c(e)R4 . 
Here c(e), C(e) are positive constants that depend only on e. 

There exists a combinatorial procedure that allows one to cut the tentacles off 
and to get an almost optimal estimate for the diameters of the modified basins. 

Theorem 1.7 (modified basins) Given e > 0, there exist open pairwise disjoint 
sets B'(a) with the following properties: 

(i) m 2 B'(a) = n; 

(ii) C = B'(a) (up to a set of measure 0); 

aeZf 

(in) m 2 (B(a) f] B'{a)) ^ tt - s; 
(iv) for any zeC, and any R ^ 2, 

P{diam(5'), > R} ^ C(e)e~ cR4 ^ R) ~ 3/2 . 
Here (B') z is the modified basin that contains the point z. 



5 



The estimate in item (iv) is not as good as the tail estimate e cR4 ( l °s R ) 1 that can 
be obtained by modification of the proof in [To! part II], but it comes fairly close. 

Now, a few words about the tools we use in the proofs. First of all, it is the "almost 
independence" of the localizations of a G.E.F. to distant disks (Theorem 13.21) . which 
may be useful in other problems as well. In the proof of the long gradient curve 
theorem, we use lower bounds for the determinants of large covariance matrices of 
some Gaussian complex random variables. These bounds are proved in Section 5. The 
proofs of the distance to the sink theorem 11.51 the diameter of the core theorem 11.61 
and the modified basins theorem 11.71 are based on a version of the length and area 
principle (Proposition 18.21) . 
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2 Preliminaries 

2.1 Basic facts about complex Gaussian random variables 

We fix some probability space (Q, P) and some (very big) family {Ej}j e j of inde- 
pendent standard complex Gaussian random variables on that probability space (i.e., 
the density of on the complex plane C is ^e - ^' ). Every complex Gaussian ran- 
dom variable in this paper will be just a (possibly infinite) linear combination of Hj 
with square summable coefficients. Such a complex Gaussian random variable £ is 
standard if E|£| 2 = 1. 

A useful remark is that if rjk are standard complex Gaussian random variables 
and ak G C satisfy J2k \ a k\ < +oo, then J2k a kVk can be represented as at] where 
^ a ^ and 77 is some standard Gaussian random variable. 

We shall start with simple probabilistic estimates. 



6 



Lemma 2.1 Let rjk be standard complex Gaussian random variables (not necessarily 
independent). Let > 0, a = J2k a k- Then, for every t > 0, 



k 



i 2 

a k \r]k\ ' 



Proof: Without loss of generality, a = 1. We have 

^ e"^ 2 ^a fc Eexp {||?7fe| 2 } = e~^ 2 Eexp {||?7| 2 } 

k 

where t] is a standard complex Gaussian random variable. But 

"+0O 



Eexp{±|r/| 2 } = i ^ e~^ 2 dm 2 {z) = 2 



re 2 r2 dr = 2 



□ 



Lemma 2.2 Lei {Ci}i^?sgrt ^ e complex Gaussian random variables, and let T = (7^) 
6e their covariance matrix; i.e., 7^ = E£j£y. Suppose det T ^ 1. T/ien 

P{ |&| ^ e, 1 < z < n} < e 2n . 

Proof: The joint density function of the variables £j is 

_ e -Cr-*M> < ^-n _ 



7r n det T 

Thus 



|fl|<£,...,|£„|<£ 



□ 



Now we want to elaborate on the well-known fact that a family {£;}ig,r of complex 
Gaussian random variables is independent if and only if the covariances E£j£j vanish 
for % ^ j. 

Lemma 2.3 Let £/. be standard complex Gaussian random variables whose covari- 
ances •jij = E£j^- satisfy 



^2 iTwl ^ °" ^ ^ for alii. 



Then = + 6^77^ where (k o,re independent standard complex Gaussian random 
variables, rjk are standard complex Gaussian random variables, and bk G [0, a). 
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Proof: Note that ||M|| = supj^j defines a norm on matrices M = (m^) (more 

i 

precisely, it is the norm of M as an operator in £°°). Now let T = (7^) be the 
covariance matrix of the family Note that T = I — A where I is the identity 
matrix and ||A|| ^ a. Then T~2 — I — A' with ||A'|| ^ a. Indeed, using the Taylor 

series (1 — z)~^ = 1 H — z + otiz 1 and observing that \ai\ < 1 for all t ^ 2, we get 

02 



lA'l 



£>2 £>2 



a 2 



1 -(7 



^ a. 



It remains to put Cfc = (r 2 ) fc j6j, hVk = A fc, 7 Cj- D 



2.2 Operators T w and shift invariance 

The main thing we need from / and U is their shift invariance. It is literally true 
that U is shift invariant (as a random process) but U is a little bit less convenient 
than / to work with because, firstly, it is not a Gaussian process and, secondly, it 
has singularities. The random function / itself is not shift invariant, but there is a 
simple transformation that makes a shift of / a G.E.F. again. 

For a function / : C — > C and a complex number tofC, define 



T w f(z)=f(w + z)e- zw e->f 



Lemma 2.4 Let /: C — » C be an arbitrary function and let w G C. Let U(z) 
t;\z\ 2 and let U w (z) = log \T w f(z) \ — || 



log \f(z)\ - \\z\ 2 and let U w {z) = log \T w f(z)\ - l\z\ 2 . Then U(w + z) = U w {z). 



Proof: 

,)=\og\T w f(z)\- 

log + z)\ - \\w + z\ 2 = U(w + z) 



U w (z) = log \T w f(z)\ - l\z\ 2 = log \f(w + z)\ - Rezw - \\w\ 2 - ±\z\ 2 



□ 



Lemma 2.5 For any w', w" G C, we have 

Tl f Hmw'w" rp rp f 

w'+w"J — e 1 w' 1 w"J ■ 



8 



Proof: 



(T w ,T w „f)(z) = (T w „f){w' + z)e-™'e->'f 

= /( y + w > + z)e -{*+ w >)^ e -h\™r e -^ e -hw\ 2 
— f( w ' + w" + z y- ilmw '^ e - z W+w 77 ) e -\W+™"\ 2 



-ilmw 



(T u; / +u ,///)(z) . 



□ 



Lemma 2.6 Lei f be a G.E.F. Then T w f is also a G.E.F. 

Proof: It suffices to check that the covariances of these two complex Gaussian pro- 
cesses are the same. Recalling that, for a G.E.F. /, we have Kf(z')f(z") = e z z " , we 
get 

E(T w f)(z')(T w f)(z") = e- z ' w -^ w e-^ 2 Ef(w + z')f(w + z") 

= e - z '™-z"' w e - M 2 e (™+z')(w+z") __ gz'i 77 __ E/(2; / )/(z // ) . 

□ 

Note that one can give another proof of this lemma using the fact that the func- 

z n 

tions — = form an orthonormal basis in the Fock-Bargmann space 7i (that is, in the 

closure of the analytic polynomials in the weighted space L^-e - ^' dm-i) ) and that 
T w is a unitary operator on 7i. 

Lemma l2l6l together with Lemma EH immediately imply that the random potential 
U corresponding to a G.E.F. / is shift invariant (as a random process). 

3 Almost independence 

Let now T w f(z) = ^^-z k be the (random) Tayl or series of T w f at 0. Lemma [231 

implies that, for a fixed w G C, Ck( w ) are independent standard Gaussian random 
variables, but, of course, the covariances between and £j(w") may be nontrivial 

for w' w". 

Lemma 3.1 

l^jK)&K)l < 5^ j+k) e->'- w "f . 
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Proof: Let w = w" — w' . Since, according to Lemma 12.51 T w iif = e llmww 'T w T w i f ', 
the random variable £k(w") equals e lImww ' \/k\ times the fc-th Taylor coefficient of 
the function ^^(w') e -zu> e -\\w\ 2 ^ jj ence ^he absolute value of the covariance 

in question is just \fk\ times the absolute value of the fc-th Taylor coefficient of the 
function 77^ (n> + z)^e~ zw e~^ w ^ . According to the Cauchy inequality, this coefficient 
does not exceed 



1 - fc 



p max 



If. \z\=p 

for any p > 0. Choosing p = j\w\, we get the estimate 



< -^= p - k (\w\ + p y e pH e ->\ 2 
Vr- 



Exchanging the roles of w' and w", we get the symmetric inequality 

k 



|E£>')&K)I < ^ (^j \w\ k -ie-* 



Taking the geometric mean of these two estimates, we get the statement of the lemma. 

□ 

Our next aim is to show that the G.E.F. T Wj f can be simultaneously approximated 
by independent G.E.F. in the disk \z\ ^ r if all distances between the points Wj are 
much greater than r. More precisely, the following statement holds. 

Theorem 3.2 (almost independence) For every N > 0, there exists A = A(N) > 
such that, for all r > 1, and for all families of points Wj G C satisfying \wi — Wj\ #S 
Ar, % ^ j, we can write 

Twjf fj ~l~ hj 

where fj are independent G.E.F. and hj are random analytic functions satisfying 
F{mwL\hj(z)\ >e~ Nr2 } < 2exp {-| exp{iVr 2 }} . 

Proof: Fix two constants A ^> B ^> 1 to be chosen later. Consider the standard com- 
plex Gaussian random variables £,k{wj) with k ^ B 2 r 2 . We want to apply Lemma [2731 
To this end, we need to estimate the sum of covariances \E£f,(wj)£,e(wi)\. 

(kj): (kJMH) 

Recall that F£k(u>j)£e(u>i) = if i = j. For j ^ i, we can use Lemma I37TI which yields 

\EUwj)Um)\ ^ (B 2 r 2 + l)5 B2r2 e->'- Wi|a . 
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It remains to estimate Y^j-j^i e ■ Let fx be the counting measure of the set 

{vjj}j^i. We have 



fx(D(wi,s)) ^ 



0, if s < Ar , 
9A~ 2 r~ 2 s 2 if s>Ar. 



(The second estimate follows from the observation that the disks D(vjj, \ Ar) are 
pairwise disjoint and contained in the disk D(wi, s + \Ar) C D(u>i, |s) if \wj — Wi\ ^ s 
and s ^ Ar). Now, write 

e ->^ 2 = jj e -il*-«*l a ^(z) = ^-^^(^K, s)) rfs 

/•oo 3 

<; 9 / A" 2 r- 2 -e-^ 2 rfs = 9 (1 + 4A-V- 2 ) e-i AV < lOe"^ 2 , 



3 3'fr 



provided that A ^> 1 and r > 1. Using this estimate, we finally get 



provided that A ^> B 1 and r > 1. Applying Lemma 12.31 we conclude that 
£,k{wj) = Ck(wj) + bkjT]k(wj) where Ck(wj) are independent standard Gaussian random 
variables, 7]k{wj) are standard Gaussian random variables, and ^ 6jy $C e~5 A2r2 . 

For k > B 2 r 2 , let Cfc( w j) be independent standard complex Gaussian random 
variables that are also independent with Ce( w i) f° r an ^ ^ -B 2 r 2 and for all z. Put 



fc>0 vfc. 



and 



hj(z) =T Wj f(z)-fj(z) = - &M~7p + £ ^ fc W7P + S ^jVkiwj 

k>B*r* ^ k>B 2 r 2 ^ fc<B 2 r 2 V 

The G.E.F. /j are, clearly, independent and all we need to do now is to show that hj 
are small in the disk \z\ ^ r. We shall use Lemma [2. II It reduces our task to that of 
estimating the sum 

2 V + V T + T -£=. 

z-^i rn jL^i k j rn /-^ rn /-^ /u 

k>B 2 r 2 V fc ^ B 2 r 2 V /V. fc>S 2 r 2 V k<B 2 r 2 V 
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Note that in the series — = the ratio of each term to the previous one equals 

k>B 2 r 2 vk\ 

<-^- = 4-<iif5>2. Hence the sum does not exceed twice the first term of 



v/fc+I ^ Br B 2 

the series, which is 



VW- " \VkoJ \ Br J ^ V B 

provided that B > e^fe (here ko is the smallest integer bigger than B 2 r 2 ). On the 
other hand, the Cauchy - Schwarz inequality yields 



k 

V ^ VB 2 r 2 + 1 

k<B 2 r 2 V 



V — = VB 2 r 2 + le* r 
^ k\ 



Thus, the sum we need to estimate does not exceed 



provided that A > B > VN. 

It remains to apply Lemma [2.11 with a ^ e~% Nr * , t = e~ Nr " . □ 

4 Size of the potential U 

First, we estimate the probability that the maximum of the random potential U over 
the disk of radius p is large positive. 

Lemma 4.1 For p ^ 1 and M > 0, 

P { max U{z) > M } ^ Cp 2 e~ ce2M . 

Proof: Since U is a stationary process, and since the disk {\z\ ^ p} can be covered 
by Cp 2 copies of the unit disk, it suffices to show that 

9 A/f 

F{m&xU(z) > M] < Ce~ ce . 
But this probability does not exceed P { max \f(z)\ > e M }, which, in its t urn, does 

|z|<l 



not exceed P { — ^=|Cfe| > eA/ } • Estimating the latter probability by Lemma [2. 11 
k vk\ 

we get the desired result. □ 
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Lemma 4.2 Suppose that p ^ 1. T/ien 

P{ max \f{z)\ ^ e" 3p2 } < e" 8p4 . 

Proof: Assume that max \f(z)\ ^ e~ 3p2 . Then by Cauchy's inequalities for the Taylor 
coefficients of analytic functions, we have 

\l n ! n n i 2 o 

|e„| ^max ^ — e" 3 ' 2 , n = 0, 1, 2, ... . 

The probabilities of these independent events do not exceed (np- 2 ) n e- 6p . Thus 
{maxl/O)! < e" 3p2 } ^ JJ [(np" 2 )" e - 6p2 ] < (2 2p ' < e~ 8p4 . 



P 



□ 

Theorem 4.3 Given (3 > 0, suppose that p is sufficiently large, and that log 2 p ^ 
M ^ p 2 . Then the probability of the event 

jthere exists a curve 7 C p 4 D with diam(7) ^ /3p such that xaaxU < — m| 
does not exceed e~ cpM3/2 with the constant c depending on (3. 

Recall that by D we denote the unit disk in the complex plane centered at the origin, 
tO is the disk of radius t concentric with D. 

Proof: We fix a sufficiently small constant a < min (|, f) and cover the disk p 4 D by 
the disks Dj = D(wj,a\/M), j E J, with bounded multiplicity of covering. Clearly, 

Suppose that there exists a curve 7 C p 4 D with diameter at least f3p and such 

that maxt/ < —M. Note that if 7 enters the disk D~, then it must exit the disk 

7 

2Dj = D(wj, 2aV / M); otherwise, Aay/M (the diameter of 2Dj) is larger than f3p (the 
diameter of the curve 7), which is impossible due to our choice of a. 

Let A be the constant corresponding to the value N = a~ 2 in the almost inde- 
pendence theorem 13.21 Having the curve 7 and the constants a and A, we choose a 
sub-collection of well-separated disks Dj, j G J* , with the following properties: 



\wi — Wj\ ^ 2AayM for j 7^ %\ 
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Figure 2: The curve 7 and the disks Dj = D(wj, ay M) 



the curve 7 enters each of the disks Dj; 



(3p 



2AaVM + 2aVM 



(3p 



2(A + l)aVM 

By \x~] we denote the least integer n ^ x. 

Applying Theorem 13.21 with r = 2a\fM, we get T Wj f = fj + hj, j E J* , where fj 
are independent G.E.F. and 

P{ max^ \hj\ > e" 4M } ^ 2exp -^exp(4M) 



2aVMl 
-AM 



If max max \hj\ ^ e , then, for z G (7 — fl 2avMD, and for big enough M, 



-M+2a 2 M 



+ e 



-AM 



< e 2 



M 



Now, we introduce the independent events We say that the event oc- 

curs if there exists a curve 7-,- that connects the circumferences = aV'M j and 

= 2aV / Mj such that < e~5 M everywhere on jj. 



Claim 4.4 // the constant a is small enough, then P { (*.,•) } ^ e 



-cM- 



Proof of Claim \4^4\ Consider the function log \fj\ subharmonic in the disk 2avMI 
By Lemma 14. 1^ throwing away an event of probability less than 



Ca 2 Me-^ a < e~ cM2 
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we have 



and hence 



max 

z£2a\fM\ 



logl/iWI 



V 



< 2a 2 M 



max log \fj\ < 4a 2 M. 

2a v / MD 



The curve jj connects the circumferences {\z\ = ay M} and {\z\ = 2ayM}. Hence 
its harmonic measure with respect to (^2ay/M3^ \ jj is bounded from below by a 

positive numerical constant Co uniformly in the disk a-y/MD (this well-known fact 
follows, for instance, from (2j Theorem 3-6]). Thus 

max log | 1 ^ 4a 2 M — —M < -yM, 

2 4 



if the constant a was chosen so small that a 2 < jqCq. Then \cqM > 3(a\^M) 2 and we 
can apply Lemma 14.21 to the function fj in the disk av^MO. The lemma yields that 
the probability that happens does not exceed e~ cM2 . □ 

We conclude that the existence of a curve 7 satisfying the assumptions of the 
theorem implies existence of a subset J* <Z J with j^J* 
at least one of the following happens: 



§p 



2(A+l)aVM 



such that 



(i) for j G J7"*, the independent events occur with 7j = (7 — Wj) PI 2av MB 



(ii) max max > e 

J'GJ"* 2a\fM D 



-4M 



In the case (i) , the probability is bounded by 



Since \ ) 



^ Cp 8 , and #J* ^ Cp, the first factor is bounded by e plosp . The second factor 
does not exceed e ~ cM3/2 p_ Since logp ^ M 1 / 2 , the whole product is bounded by 
e" cM3/2p . In the case (ii), the probability does not exceed 



#J* ■ (*^\ ■ 2e" ceCM < e c ^"-^ U < e— cM 
which is much less than e - c P M ' i/2 . This completes the proof. □ 

5 Determinants of covariance matrices 

In this section, we estimate from below the determinant of the covariance matrix of 
the complex Gaussian random variables {f'(zi) — Zif(zi)}i^i^ n . This estimate will be 
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used in the next section when we apply Lemma 12.21 to the proof of the long gradient 
curve theorem 1 1.41 

To warm up, first, we estimate the determinant of the covariance matrix of random 
variables which has a simpler structure. 

Lemma 5.1 (the 1st determinant estimate) Let {zi}i^i^ n C C be a well- separated 
sequence; i.e., for some A > 0, 

\zi — Zj\ ^ A|z — j\, l^i,j^n, 



and let T = (7^) where 7^ = Kf(zi)f(zj) = e ZiZj , 1 ^ i,j ^ n. Then 

detr> (cAv^) re(re_1) • 

Proof: Without loss of generality, we suppose that n ^ 2 (if n = 1, the statement is 
obvious). Since 

00 z k z k 



E 

k=0 



Id V¥. 



we have T = AA* with the matrix 



.4 



V vT< V». 



4 


z l 






z 2 


z k 







Hence, by the Cauchy-Binet formula, 

detr = \ m M)\' 



where the sum is taken over all principal minors m t (A) of the matrix A. We use only 
one principal minor 



m (A) 



1 -43= 
1 Vv. 



1 z n 

1 Vv. 



'l!2!...(n-l)!n ( 



Zj Z{J . 



Since the points z±, z n are well-separated, we get 

\m (A)\ ^ A^Vl^!...^-!)! 
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Since k\ ^ k k e k , k ^ 1, we have 

(71-1 
^(fclogfc-fc) 
k=l 



n(n — 1) 

^ exp I / x log xdx — n log n 



^ exp n ~ 11,2 ~ n l°& n 

and 

\m (A)\ 2 > (Av^) n(n - 1} ■ e - n2 ~ n ^ n > (cX^f^ , 
completing the proof of the lemma. □ 
In the second estimate, we fix the parameters n G N and r = By/n where B ^> 1. 

Lemma 5.2 (the 2nd determinant estimate) Let {^ji^^n C C be a collection 

r ^ ^ 

of points such that \Zi — zA ^ - i-j, and Zet min ^ r. Lei T &e i/ie covariance 

n i 

matrix of the complex Gaussian random variables & = /'(Zj) — z~if(zi), 1 ^ i ^ n. If 
B is sufficiently big, then det T ^ 1 . 



The idea of the proof of this lemma is similar to that of Lemma I5TT1 though the proof 
is more involved due to a more complicated structure of the covariance matrix. 

Proof: First, we compute the values 7^ = E^^: 
Claim 5.3 7^ = (1 — |z, — Zj\ 2 ) e ZiZj . 
Proof of Claim 15.31 : 



E(f'(z) - zf(z))(f'(w) - w f(w)) = (d z -z){c\- w) Ef(z)f(w) 

= (d z - z) (du, - w) e zw — (1 + wz - zz - ww + wz) e zw 

= (1 - |^ - w| 2 ) e zW 



□ 



Now, we suppose that n ^ 2 (if n = 1, the statement is obvious) and factor the 
matrix V. We have 



la = (d Zi - Zi) (d 8 . - z 3 ) 



00 zf z k 



v k\ Vk\ 



£ 

k=0 
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Put 



.4 



■zi 



V 





(2- 




(fc- 




yffl 










l-\z n \ 2 


(2- 


| -Sri | 


(fc- 




yffi 




\/2! 




sfk\ 



■J 



Then T = AA*, and by the Cauchy-Binet formula, 

n+l 



t=i 



where the sum is taken over n + l principal minors M t (A) of the matrix A: 

{k + t _\ Zi \2 )z k + t-l 



M t = det 



t= 1,2, ...n + l. 



To estimate the sum of the squares of these determinants, we introduce the determi- 
nants of simpler structure: 



fit = det 



(k + t- \zi\ 2 )z, 



i 



t = 1,2, ...n + l 



l<i,fc<n 



Claim 5.4 For 1 < t < n + 1, |M t | ^ 



Proof of Claim \5.4[ follows by a straightforward estimate of the ratio 



{\Z\\ ■ \z 2 \ ■ ... ■ \Zn\f ■ 



1! 



2! 



ni 



Since \zi\ ^ r, we have {\z\\ ■ \z 2 \ 
n+l, 



(t + 1)! (t + 2)\ "' (i + n)!" 
^n|) ^ r * n - F° r each integer i between 1 and 



1! 



2! 



??.! 



(t + l)\ (t + 2)\ "' (t + n)\ " V(* + ™) ! 



n/2 



(n + 1) ... (n + t)J \n + t 



n/2 



nt/2 



>'s) 



nt/2 



Thus 



if 

3rt 



3 



^ 1, 



provided that B ^ \/3. This proves the claim. □ 

Now, we complete the proof of Lemma 15.21 Observe that m is a polynomial of 
degree n int. We use a version of the pigeonhole principle: 
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Claim 5.5 Let P(t) be a polynomial of degree n with the leading coefficient a. Then 
max{|P(t)|: t e {1,2, ...,n+ 1}} > \a\2~ n . 

Proof of Claim 15.51 ' We have 



max |P(i)| = |a| • max \t — T\\ ... \t — r, 

ie{l,2, ... ,n+l} <e{l,2, ... ,n+l} 



n I ) 



where ri, r n are the zeroes of P(t). We have n+1 disjoint —neighbourhoods of 
the points 1, 2, n + 1 in C. At least one of them is free of the zeroes of P(t). 
Hence at the center of this neighbourhood, the absolute value of P cannot be smaller 
than |a|2 _n ; whence the claim. □ 

We apply this claim to the polynomial \i t . Its leading coefficient equals 



lim t n /if = det 



'k\J 1 



We've already estimated this determinant in the proof of the model lemma 15.11 We 
get \a n \ > (cAv^) n(n " 1)/2 with A = - and n = (r/B) 2 . Then \a n \ ^ (cP)^"- 1 )/ 2 , and 



x w(n-l)/2 



max \u t \ ^ \a n 2 n ^ -cB 
te{i,2, ...,7i+i} \2 

whence 

n+1 n+1 ,^ \ 71(71— 1) 

detr>J]|M t | 2 ^^|/i t | 2 > (-cB) 

t=i t=i / 

provided that the constant P is chosen sufficiently large. This completes the proof of 
Lemma 15. 2[ □ 



6 The long gradient curve theorem 

Till the end of the proof, we fix 5 = ^. Everywhere below we shall assume that 
R 3> 1. In the proof, we work with three scales: starting with the macroscopic R- 
scale, we move to the intermediate y/\ogR- scale, and then to the microscopic -^-scale. 

6.1 Bad squares 

Suppose that there exists a gradient curve T connecting 8Q(0, R) and dQ(0, 2R). Due 
to Lemma \A. 11 we can assume that U ^ R s everywhere on Q(0, 2P), hence on T: the 

probability of the opposite event does not exceed CR 2 e~ ce . Suppose that there is 
a point on T where U = —R s . Since T is a gradient curve, if such a point exists, then 
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it is unique. This point splits T into two parts: I\ where U < —R s , and T 2 where 
U ^ —R s . If U ^ — i? 5 on T, then one of these parts is empty. One of the curves T 1 , 
r 2 must connect either dQ(0,R) with dQ{Q,y/2R), or dQ(0, V2R) with dQ(0,2R). 
If this is the curve Fx, then its diameter is larger than cR. By Theorem 14.31 the 

probability of this event does not exceed e~ cR 3 , and we are done. 

Thus the proof boils down to the case when the gradient curve T connects dQ(0, R) 
with dQ(0, V2R) and — R s ^ U ^ R s everywhere on T. In this case, 

\VU(z)\ \dz\ ^ 2R S ; (6.1) 

that is, the gradient VU is small in the mean on T. We will not use anymore that T is 
a gradient curve; starting this moment, it is an arbitrary curve connecting dQ(0, R) 
with dQ(0, y/2R) such that (16. ID happens. 

We take r = —y/S log R and fix the standard partition of the complex plane C 

into squares Q(wj,r) with side length 2r. Let J be the set of indices j for which the 
square Q(wj,2r) is entirely contained in the "square annulus" Q(0, y/2R) \ Q(0,R). 
Note that #J ^ (R/r) 2 . 

Definition 6.2 (bad squares) Let j G J. We shall call the standard square Q(wj,r) 
bad if there exists a curve jj joining dQ(vjj,r) with dQ(vjj,2r) such that 

j \VU(z)\ \dz\ < rR 25 - 1 . (6.3) 

We shall call the square Q(vjj, r) good if it is not bad. 

By T C J we denote the family of all indices j such that the square Q(vjj,r) 
intersects the curve T. 

Lemma 6.4 At most 8i? 1_<5 jr of the squares {Q(wj, r )} JgJC - are good. 

Proof: Let N denote the number of good squares Q(vjj,r). By we denote a 
connected part of T fl (Q(wj,2r) \ Q(vjj,r)) that joins dQ(wj,r) with dQ(wj,2r). 
Since almost every point of the curve T belongs to at most 4 squares Q(wj,2r), we 
can write 

NrR 25 - 1 < J l Vf/ (^)l \ dz \ < 4 ^ l Vf/ WI M 2 I < 8i?<5 , 

j: Q(wj,r) is good ^ J 

whence the estimate. □ 

The immediate consequence of Lemma 16.41 is that the existence of a curve T 
connecting dQ(0, R) with dQ(0, \/2R) such that (16. lft happens implies the existence 
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cR 

of a family T of squares Q(wj, r) of cardinality L ^ — and a subfamily T' C of 



bad squares of cardinality at least L 



8R 



1-8 



^ -L in that family. 




Figure 3: The gradient curve T generates separated bad squares 

Let A be the constant corresponding to N = 4<5 _1 in the almost independence 
theorem 13 .21 Let J' C J satisfy \wi — Wj\ ^ 4Ar for i,j G J\ i ^ j. According to 
Theorem 13.21 applied to 4r instead of r, we can represent T w .f as fj + hj where fj are 
independent G.E.F. and all the functions hj are small in the disk \z\ ^ 4r. We set 

fL = {max max \hAz)\ > R~ } . 

jej' |*|<4r J 

Then for any j G J', P{ max > R~ 4 } < 2e~^ 4 (recall that e " 167Vr2 = iT 4 

for our choice of r and N). Therefore, P{f2* } ^ 2R 2 e~^ Ri < e~ cRi . 

The next proposition is the central part in the proof of the long gradient curve 
theorem: 

Proposition 6.5 There exist events Qj with F {Qj } «C e~ cr4 depending only on fj 
(and, thereby, independent), and such that, for any j G J' , 

{ Q(wj, r) is bad } C Qj U f2* . 
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Now, using this proposition, we complete the proof of the long gradient curve 
theorem. We choose a family T" C T' of cA~ 2 r~ l R "4v4r-separated squares" (that 
is, all the distances between the centers of these squares are not less than 4Ar), and 
discard the rest of T' '. From Proposition 16.51 we see that the probability that a given 
subfamily T" C J of cA~ 2 r~ l R squares is bad does not exceed 

provided that R 3> 1. At last, we have at most 

{#jf T " <C (CR 2 ) CR <e CRh * R 

ways to choose T" in J . This does not harm the previous upper bound. Hence the 
long gradient curve theorem is proved (modulo the proposition). □ 

6.2 Proof of the proposition 

Assume that the event does not occur. Then T w .f = fj+hj where max \ hA ^ R~ . 

3 J J 4rD 

We fix j G J' and aim at building an event Qj depending only on fj of probability 
P { Qj } ^ e~ cr and such that, if the square Q{vjj, r) is bad and hj is small as above, 
then Qj must occur. To simplify the notation, we set w = Wj. 

Fix the partition of the complex plane C into standard squares with side length 

2_ 
R' 

Definition 6.6 (black squares) We shall call a standard square Q((, j^) C Q(w, r) 
black if inf |W| sC R 35 " 1 . Otherwise, the square Q((, 4) is called white. 

First, we check that if the square Q(C, j=>) is black (i.e., the gradient V£7 is small 
somewhere in this square), and the functions fj and hj are not too large, then the 
function fj(z) — zfj(z) must be small at the center ( — w of the shifted square. 

Lemma 6.7 Suppose that 

(i) the square Q((, j^) is black; 

(ii) max \ hA ^ R~ ; 

4rB 1 Jl 

(iii) max | £ Kit! 5 . 

4rH 

Then 

|/j(C -w)- (C^)fj(C - w)\ < R 65 ' 1 . (6.8) 
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(T w f)'(z) 



Proof: We have U{w+z) = log \T w f(z)\ — ^\z\ 2 , whence \VU(w+z)\ 
Thereby, 

Iffa) - zfj(z)\ ^ {h'jiz) - zh,{z)\ + \VU(w + z)\ {\f s {z)\ + Mz)]) . (6.9) 

Since the square Q((, j^) is black, there exists a point z such that w + z £ Q(C, and 
| V£/ (u>+z) | ^ R 36 ^ 1 . The other terms on the RHS of (16. 9p are readily estimated using 
assumptions (ii) and (hi) and Cauchy's inequality for the derivative of an analytic 
function. We get 

\f 3 (z) - zfj(z)\ ^ 2riT 4 + R 35 - 1 (R 5 + iT 4 ) < R 55 ' 1 . 

It remains to replace z by Q — w on the LHS. 
By Cauchy's inequalities, 



max I f'l < r~ 1 R d , and max I f'\ 2r- 2 /T . 

Q(0,2r)' 51 Q(0,2r) ' 

Hence the operator norm of the differential of /j(^) — zfj(z) does not exceed 

2r~ 2 rR 5 + 2rr~ 1 R 5 + 2^ < 5i2 5 
everywhere in Q(0, 2r). Since \z — (£ — u;)| ^ and i? ^> 1, we are done. 



□ 




Figure 4: The curve 7 and a sequence of black squares it generates 

Assume that the square Q(w, r) is bad; i.e., there exists a curve 7 joining dQ(w, r) 
with dQ(w, 2r) such that 



|VZ7(*)| < ri? : 



2<5—l 
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We fix an integer n = (r/B) 2 with B 1. For any t G [0, 1], we put 

(i-l)+t 
rAt) = r H r , 1 ^ z ^ n . 

For each t, the squares dQ(w, ri(t)) form a "chain of n fences" , and the curve 7 crosses 
this chain at least n times. It may happen that, for some value t, the gradient V£7 is 
not small at most of the crossing points, or even at all of them. However, as we shall 
see, for a large subset of t G [0, 1], the gradient VU is sufficiently small at n crossing 
points 7 fl dQ(w,ri(t)), 1 ^ i ^ n, to guarantee that the corresponding -^-squares 
containing these points are black. For each t G [0, 1], we denote by B(t) C {1, 2, ... n} 
the subset of those z's that at least one point from the set 7 fl dQ(w, is covered 
by a black square. By mi we denote the one-dimensional Lebesgue measure. 

Lemma 6.10 Suppose that the square Q(w,r) is bad. Then 

mi{te[0,l]:#S(t)=n}>i. (6.11) 

Proof: Let L be the measure of the set of p G [r, 2r] such that the intersection 
7 fl dQ(w,p) is contained in white squares. Then 



f \VU{z)\ \dz\ > ti?*- x L. 



Since the square Q(w,r) is bad, the LHS does not exceed rR 2S 1 , and we see that 
L < rR~ s . On the other hand, 



(in- #B(t)) dt^-- mAt G [0, 1] : #B(t) < n - l) 
Jo n 1 



2 

n r 



whence mi{t G [0,1]: #5(t) ^ n - l} ^ ^ < 1 if R > 1. Hence the 

lemma. □ 

For each t G [0, 1], consider the collection 3(£) of "configurations" 3 = {zi, ... z n } 
of n points such that each point Zi is a center of a standard square Q (zi, ^) from 
our partition that has a non-void intersection with dQ(w, ri(t)). Let us introduce the 
events 

TV( 3 ) = ( max Iffa - w) - (i^T)/-^ - w)| < i? 6 ^ 1 ) and fi^i) = \J T^), 

L <<n J J63(t) 

and estimate their probabilities. Our estimate is based on the lower bound for 
the determinant of the covariance matrix of complex Gaussian random variables 
{fj(zi) - Zifj(zi)} given in Lemma E21 
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Lemma 6.12 Given t £ [0, 1] , P{ fy(t) } < e~ cr \ 

Proof: First, we estimate the probability of the event 

Claim 6.13 For any configuration 3 £ 3(t), P{ ^(3) } < fl-2(i-W)» 

Proof: This is a straightforward combination of Lemmas 12.21 and 15.21 □ 

Next, we estimate the cardinality of the collection 3(£) (recall that t £ [0, 1] is 
fixed) . 

Claim 6.14 #3(t) < C n (Rr) n . 

Proof: For each i, there are at most CRr standard -^-squares that intersect dQ(w, 
Therefore, there are at most CRr choices for the centers z% of these squares and the 
number of the corresponding configurations 3 cannot exceed (CrR) n . Hence the claim. 

□ 

Using Claims IB. 141 and [6.131 we get 

P { Qj(t) } <: C n {Rr) n ■ iT^ 1 - 65 )" ^ C n R -{l-135)n ; 

13 r 2 

provided that R is sufficiently big. Recalling that 135 = — < 1, n — — , and 
r = — 51og.R, we see that P{fij(t)} ^ e~ cr . This proves Lemma [6. 121 □ 



Define the events 



fij = |cj £ fi: miO G [0,1]: w G fij(t)} ^ ^} 



and 



fi,- = fi' U | max I/,- 1 ^R s \. 

J J I 4rB J ) 

Note that the event Qj depends only on fj. By Lemmas 16.71 and I6.10[ 
\Q(w, r) is bad ) C fL U { max \hA ^ i? -4 ) . 

4rH> J 

Lemma 6.15 TTie probability of the event j max ^ i? 5 } does no£ exceed e~ cri . 



Proof: If ^ i? 5 somewhere in the disk 4rO, then (at the same point) the corre- 
sponding potential U is not less than 5 log/? — |(4r) 2 = itflog-R (due to the choice 
of r). By Lemma I4~T1 the probability of this event does not exceed Cr 2 e~ cRS <C e~ cri 
ifi?>l. □ 



Hence to complete the proof of Proposition l6.51 we need to estimate the probability 



of the event fi'-. 
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Lemma 6.16 P{^} ^ e~ cr . 

Proof: Define the random set A = {t G [0, 1] : Qj(t) occurs}, and let X = mi(A). 
Then, by Chebyshev's inequality, 

P(ftJ) = F(X ^ 1/2) < 2E(X) = 2^ F(£lj(t))dt < 2 max P(fij(*)). 

By Lemma [6.121 the maximum on the right-hand side does not exceed e _cr4 . □ 
This completes the proof of Proposition 16.51 



7 Proof of the partition theorem 



Set A(z) = U xx (z)U yy (z) -U 2 xy (z). 
Lemma 7.1 For z eC\ Z f , A(z) = 1 - 



f)'« 



Proof: is a straightforward computation. Since d x = d z + d z , d y = i(d z — d z ), we 
have 



dxx &zz 2d zz ~\~ dzzi &yy d zz -\- 1d zz d zz , d xy i{dzz dzz) • 
Whence 

A(z) = (U zz + 2U Z - Z + U- zz ) (-U zz + 2U Z - Z - U- zz ) + (U zz - U- zz f = 4 {U 2 Z - Z - U zz ■ U- zz ) 



l ff'Y l ff'Y l 

Taking into account that U zz — - — , U zz — - — , and C/ Z g = — , we get the 
result. □ 
Denote by Crit U the set {z G C: VC/(^) = 0} of critical points of the potential 

U. 

Lemma 7.2 Almost surely, the following hold: 

(i) each critical point of U is non- degenerate; i.e., A(w) ^ for w G Crit({7); 

(ii) the critical set Crit U has no finite accumulation points. 

Proof: Note that the probability that w = is a critical point of U is 0. At the 
critical points w ^ of U, we have /(w) = w" 1 f'{w). Hence 



AM = 1 



f" 

w—{w) - w 



w G Crit U . 
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Now, let us set f(z) = £ + £i£ + h(z), where h(z) is a random entire function 
determined by £2, £3, • • • • Then, on Crit £7, the determinant A(w) coincides with 




2 



A 1 (w) := 1 



h'(w) + £i 



Observe that, for each w with \w\ ^ 1 and each £2, £3, • • • , the set of £1 where the last 
expression is has zero measure. Thus, using the Fubini theorem, we conclude that 
for almost all £i,£2, • • • , the set {w £ C : A 1 (w) = 0} has zero measure. Let now 
g(z) = £]Z + h(z). If g is fixed (i.e., £1, £ 2 , are fixed) and w 7^ is a critical point 
of U, then £ is determined by equation 



The right hand side defines a real-analytic mapping of the punctured plane C \ {0} 
and, therefore, it maps sets of zero area in the u>-plane to sets of zero area in the 
£o-plane. Hence, for almost every choice of the independent coefficients £1, £2, the 
set of £0 for which there exist degenerate critical points of U has measure zero. Using 
Fubini's theorem once more, we get the conclusion of statement (i) of the lemma. 

Statement (ii) follows from (i). The planar map given by V£7 is real-analytic 
outside the set where U equals — 00. Note that, unless / identically equals (which 

is an event of zero probability), the gradient V£/(z) = yf{ z )) ~ z tends to 00 at 

every zero of / and, therefore, no point of U~ l (— 00) = Zf can be an accumulation 
point of Crit U. Thus, if the set Crit U has a finite accumulation point, then this 
point itself belongs to Crit U and, by the inverse function theorem, the map given by 
V£7 is degenerate at this point. □ 

It is worth mentioning that there is another way to prove statement (ii) of Lemma l?T2l 
elaborating on the fact that, if g is an analytic function and the solutions of the 
equation g(z) = z have a finite accumulation point, then g must be a Mobius trans- 
formation. 

Lemma 7.3 Almost surely, the following hold: 

(i) each oriented curve V has a starting point s(T) £ Crit([/) and a terminating point 



(ii) at any limiting point, the oriented gradient curve T is tangent to a straight line 
passing through that point. 

Proof: We refer the reader to [101 Chapter 4] for the facts from the standard ODE 
theory we use. 



9(w) ■ 



w 



t(T) £ U~ 



{-00} U Crit(*7); 
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(i) It follows from the long gradient curve theorem that, almost surely, gradient curves 
cannot escape to or come from infinity. Now it remains to observe that the limiting 
set L of any gradient curve T is contained in the set of singular points of the gradient 
flow; that is, in the set Crit U U U~ l { — oo}. Hence, by Lemma 17.21 C consists of 
isolated points. 

(ii) The critical points of U are either local maxima or saddle points. By Lemma 17.21 
almost surely, all of them are non-degenerate. The rest follows from the standard 
ODE theory: the behaviour of the integral curves in a neighbourhood of these points 
is the same as the behaviour of the integral curves for the linear ODE obtained by 
discarding the non-linear terms in the Taylor expansion of VC/. □ 

Lemma 7.4 Each gradient curve is real analytic everywhere except at the limiting 
points. 

Proof: VC/ is real analytic everywhere except on the set where U — — oo. Hence, by 
the Cauchy theorem, the gradient curves are real analytic at all points where VC/ 7^ 0. 

□ 

Now we are ready to prove the partition theorem 11.21 By the long gradient curve 
theorem 11.41 almost surely, all the basins are bounded. We call a gradient curve V 
singular if t(T) e Crit U. Note that, almost surely, every point that is not in one of the 
basins must lie on a singular curve. Moreover, with probability 1, for every compact 
K on the complex plane, there exists another compact K such that all gradient curves 
intersecting K are contained in K. (Otherwise, there exists an iV e N such that, for 
any integer M > N, there is a gradient curve connecting dQ(0,N) and dQ(0,M). 
The probability of this event is 0.) Also, a gradient curve cannot terminate at a 
local maximum of U and each saddle point of U serves as a terminating point for 2 
singular curves. This allows us to conclude that, almost surely, we may have only 
finitely many singular curves intersecting any compact subset of C. In particular, 
almost surely, each basin B(a) is bounded by finitely many singular curves and their 
limiting points, which is enough to justify the area computation in the introduction. 

□ 



8 The upper bounds in theorem 11.5 



First, we prove a useful "length and area estimate" of deterministic nature valid for 
Liouville vector fields; that is, the fields with constant divergence. Then we derive 
the upper bounds for the probability that a given point z is far from its sink a z . 
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8.1 The length and area estimate 



Consider the disk D = {\z — a\ < e\. Since VU(z) = -. + 0(1), as z — > a, 

\z — a\ z 

we can fix a sufficiently small e > such that each gradient curve hits the boundary 
circumference {\z — a\ = e} only once. This gives us a one-to-one correspondence 
between the points of the circumference {a + ee* 6 } and the gradient curves in B(a); 
i.e., the gradient curves are parameterized by the "angular coordinate" 0. 

By D(t) we denote the pre-image of D under the gradient flow of VU for time t; 
dZ 

i.e., if — = -VU(Z{t)), then D(t) = {z = Z(0) : Z(t) G D}. By A{t) we denote 

(JjV 

the area of B(a) \ D(t). Since div(V£7) = —2 on B(a) \ {a}, the evolution of the 
dA 

area is very simple: — = — 2A. This is Liouville's theorem (which follows from the 

at 

divergence theorem), see, for instance, [H §16]. 

We will need an "infinitesimal version" of this equation. The boundary dB(a) 
contains finitely many saddle points of U . By a± < ... < a s < a s+ \ = ot\ + 2ir we 
denote the angular coordinates of the gradient curves that connect the saddle points 
on dB(a) with the sink a. Take any 9 different from a±, ... ,a s , say a\ < 9 < ati+x, 
and choose 9i and 9 2 such that a\ < 9 X < 9 < 9 2 < ati+i- The gradient curves 
r(#!), r(# 2 ) must terminate at the same local maximum. They bound a "diangle" 
Y(9i,92) with the vertices at a and at a local maximum. Consider the "triangle" 
T(t; 0i, 2 ) = Y(9 U 2 ) \ D(t) and its area A(t; 9 U 2 ) = m 2 T(t; 9 U 2 ). 




Figure 5: "Triangle" T(t;0 1; 2 ) 

By Green's theorem, 

-2A(t;9 u 9 2 )= [[ AUdm 2 = [ (VU,n)\dz\ (8.1) 

where n is the unit normal directed outward the triangle. The boundary dT(t; 1; 2 ) 

Oil 

consists of parts of the gradient curves r(0 x ) and T(0 2 ), where = 0, and of the 
part / = /(i;0i,0 2 ) of the curve dD(t). If {(x(t, 0), y(t, 0)) : 9 1 ^ 9 < 2 } is the 
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equation of the arc J, then at the point 9 the unit normal n is given by 



V x e + vl 

1 f 6 ' 2 

Hence A(t; 9\, 9 2 ) = — / (U x yg — U y xo) d9, and we conclude that the area A(t; 9\, 9 2 ) 

2 Je 1 

has a smooth angular density S(t,9) = | (U x ye — U y xg). By Liouville's theorem, 
dA(t' 9\ 9 2 ) 

= — 2A(t; 9i, 9 2 ). Therefore, the density S(t,9) satisfies the same differ- 

• dS(t,9) nri , „, 
ential equation — — — = — 2b (t, 9). 

Now, we re-parameterize the gradient curve T(9) by its length / starting at the 
sink a. We treat the restrictions of the density S and of the gradient V£7 to T(9) 
as functions of the length Z; i.e., S(l) = S(t(l),9), and similarly for S7U. Note that 

— = IVC/I. We arrive at the ordinary differential equation for the density S: 
dt 

Solving this equation, we get 

dl 



S(l) = S(Iq) exp 



Denote by Z\ the point on the gradient curve T(9) that cuts the arc of length / 
from that curve. By the Cauchy-Schwartz inequality, 



|V(7| ' k V7fe ' ' / U(z,)-U(z h ) 

We arrive at the crucial 

Proposition 8.2 In the same notation as above, 

2(1 - l ) 2 



S(l) ^ S(l ) exp 



U( Zl )-U(z l0 ) 



8.2 Distance to the sink (the upper bound) 

Fix 5 G (0,2]. We define the tentacles Tr(o) of the basin B(a) as follows. Given 9, 
we move along the gradient curve T(9) in the direction of growth of the potential U, 
starting at the sink a, till we hit the point where U = —R 5 . After that, we keep on 
moving along T(9) the distance R (measured along T(9)), and then stop. The rest of 
the curve is called the ^-tentacle. The tentacles Tr(o) are the union of all ^-tentacles. 
Of course, it may happen that the tentacles Tn(a) are empty. 
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Now, we are ready to estimate the probability that \z — a z \ is large. By translation 
invariance, this probability does not depend on the choice of z, so we choose z — 0. 
Suppose that |ao| > 2R. We know that at least one of the following happens: 

(i) either the distance from to the curve T n {U < —R s } measured along T is 
less than R; 

(ii) or G T R (a ). 

(Recall that r is the gradient curve that passes through the origin.) 

In the first case, the curve 7 = r fl {U < —R s } connects the circumferences 
[\z\ = i?} and {\z\ = 2i?}. By Theorem 14.31 the probability of this event does not 

exceed Ce~ cR +1 . 

Now, we estimate the probability of the event (ii) . By translation invariance, 

7rF{0ET R (a )} = [[ P{ w G T R (a w ) } dm 2 {w) 
J Jo 

= / m 2 {w G D: w G T R (a w )}dF. (8.3) 
Jn 

Thus, we need to estimate the area of the random set {w G D: w G T R (a w )}; that is, 
the of the union of all possible tentacles within D. 

We throw away three exceptional events. Let Qi be the event that there exists a 
gradient curve connecting the circumferences {\z\ = 1} and {\z\ = R 2 }. By the long 
gradient curve theorem, Pj^! } ^ e~ cR . If Qi does not occur, then \a\ < R 2 , for 
any basin B(a) that intersects the unit disk. Let f2 2 be the event that there exists 
a gradient curve connecting the circumferences {\z\ = R 2 } and {\z\ = 2R 2 }. Again, 
P{fi 2 } < e~ cR . If fii and Q 2 do not occur, then any basin B that intersects the 
unit disk D is contained in the disk 2R 2 B>. Recalling that each basin has area n and 
comparing the areas, we see that the number of such basins does not exceed 4R 4 . 
At last, we exclude the event f2 3 = {max U > R 6 }. By Lemma EU P{fi 3 } < 

CR 4 e~ ceR < e~ cRi if R is big enough. 

Now, after throwing away these three events, we can estimate the area of the 
random set [w G D: w G T R (a w )}. First, we bound the area of one tentacle T R (a). 
Since U ^ R s everywhere in B(a), for each ^-tentacle, we can apply the length 
and area estimate from Proposition 18.21 with I — l ^ R and U{z{) — U(zi ) ^ 2R 5 . 
Integrating over 9, we get 

m 2 l R {a) ^ m 2 -D(a)e = Tie 

The number of tentacles coming from different basins and hitting the unit disk D does 
not exceed 4i? 4 . We conclude that if the events 1 ^ i ^ 3, do not occur, then 
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m 2 {w EB>: w E T R (a w )} ^ AirR^e r2 6 . In view of (18.31) . we see that the probability 
of the event {0 E Tr(clo)} is bounded by e~ cR2 if R ^> 1. 
Thus, 

P { |ao| > 2R } < P { diam (r n {U < -R 8 }) > R } + P { E T R (a ) } 

< Ce~ cRl+¥ + Ce- cR2 " S . 
Choosing S = |, we complete the proof. □ 



9 The lower bounds in theorems 11.31 and 11.5 



The proofs of the lower bounds for the diameter of the basin and the distance to the 

z n 

sink are based on the same idea. The function —= has a singular line l\z\ = \/n\ 

\/n\ 



where the gradient of its potential vanishes. Then after any analytic perturbations 
small in the annulus { | z — y/n\ ^ 2 } , this annulus still contains plenty of long gradient 
curves and of points that are far from their sinks. 

9.1 Diameter of the basin (the lower bound) 

We choose a big R ^> 1 such that n = R 2 is an integer and consider the function 

z n ( z \ 

F(z) = —== ( 1 H ) in the domain 

v^! V lORJ 

[ 7T. 1 

D=<{zeC:i2-l<|z|<i2+l,|argz--|< — 
Note that, for the corresponding potential U, we have 



b' n 
VU{z) = —{z)-z="-z + 



F 2 MXi + lk 

Since the vector | — z is purely radial and the sine of the angle between the vectors 

I + and z is at least | for z 6 O, we see that the angular component of — V£7 is 
oriented counter-clockwise and its size is at least in T>. Also, the gradient field 
—VU is directed outside the domain D on the boundary arcs {|z| =R±1, | argz — 

I I < jq}, and the radial component of — VU is at least 1 on these arcs. 

Thus, there is a gradient curve that starts at the right boundary interval { arg z — 
| = — ^j, R — 1 < \z\ < R + 1} , and hits the point iR. Thereby, its diameter must 
be at least ^. 

This conclusion will be preserved if, instead of the function F, we consider its 
analytic perturbation F + H with H satisfying |^| ^ R~ 2 in the annulus {i? — 2 ^ 
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Figure 6: The field —VU in T> 



^ R + 2}. Indeed, the absolute value of the perturbation of VU the function H 



l-R- 



2 < for R 3> 1, which is too 



creates in V is only t—L _ _ = ^ +( ' H/ ' F) < 

small to change anything in the above picture. 

Now it remains to estimate from below the probability of the event that a G.E.F 
/ is such a perturbation of F. 

f 



Lemma 9.1 IfR ^> 1, then P{ max 

H-2^|z|<R+2 

Proof: We write f(z) = F(z) + H(z) where 



1 



< _ 1 > p -CR{\ogRf' 2 



H(z) 



E 

k : k^=n,n+l 



+ (6. - 1; 



+1 



10R 



y n+l 



n 



1 ! 



1 lzl n 

Since ^ --LL 



in the annulus R — 2 ^ |z| ^ i? + 2, it is enough to estimate 



from below the probability of the event that 



max 

R-2^|z|^R+2 



E l&i 

fc: fc^n,n+l 



I k—n 



+ |£n-l| 



n+1 



10/2 



Vn + l 



< R~ 



say Now, let us handle £ n and £ n+ i first. We just demand that the corresponding 
terms be both less than R~ A . It is not hard to see that the probability of this event is 
about R~ 16 . We may neglect it since the factor R~ 16 does not affect the lower bound 
g-ORQogij) 3 / 2 we are ^ r yj n g ^ g e ^_ ^he remaining sum can be estimated as 



E 

k : k<n 



M^=(R-2f-+ E l^l^ + 2 ) 



\k—n 



k: k>n+l 



We shall show how to estimate from below the probability that the second sum is less 
than R~ A . The estimate for the first sum is very similar and we omit it (note that 
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the corresponding events depend on different and, therefore, are independent, so 
the probability that both sums are small is just the product of the probabilities that 
each of them is small). Let k = n + m, m = 2, 3, . . . . We choose some big constant 
A ^> 1 and split the sum into two: Si = and S2 = . We shall 

2<m^AR v / Iog~H m>AR y / TogR, 

show that the probability that IS2I < R~ 5 is very close to 1 and the probability that 
I & I < R^y^gR is at least e ~CR(io g R)^\ 

To estimate S2, we would like to use Lemma I2TT1 . To this end, we need to estimate 
the sum 

(R + 2) m 

^ J(n + Wn + 2) ...(n + m) ' 

m>AR^/\ogR V V A / V I 



Note that, starting with m = n, the terms in this sum decay like a geometric progres- 
sion, more precisely, the ratio of each term to the previous one is < ^= = 

^7=^ < I if R is large enough. Thus, it is enough to estimate the sum over m such 

that ARy/\og R < m «C n. Now, for k = 1, 2, . . . , n, we have n + k ^ ne^ . Thus, the 
m-th term of our sum does not exceed 

9 \ m ™ „ 2 , „2 _2 

_ \ T r A: 2m m 2 m 



R 



k_ 2m _ m 2m m _ m 

e 4 « < e R e »™ = e « e zr 1 ^ Ce mh 7 , 1 ^ m ^ n , (9.2) 



fc=i 

and the whole sum does not exceed 



2 t 2 



E_ m / _ t 

e i6A 2 <C C / e i6« 2 (it 

m>AR\/log R 

f 00 22 2 

sC CR / e~ie ^ CRe~^ losR = CR 1 ^^ < iT 6 

if i? is large enough. Thus, according to Lemma I2~T1 the probability that IS2I < R~ 5 
is very close to 1 and, at least, greater than |. 

As to Si, we just demand that each term in Si be less than R~ 6 (then |5i| < 

(R+2)m E3 

AR~ b J\og R). Since the coefficients . — 1 ; = ^ C, it is enough to demand 

V ^ ' x /(n+l)(n+2)...(n+m) to 

that |£ fc | < C" 1 ^ 6 forn + 2^/;;^n + ARy/log R. But the probability of this event 
is at least ( c r-^)Arvwr ^ e -CR(\ og Rf/\ Thig proveg the lemma n 

Thus, with probability e - CR ( l °s R ) 3/2 ; the point z = iR belongs to a basin of diam- 
eter greater than y|. It remains to note that, due to shift invariance of U, the same 
is true for any other point z on the complex plane. This proves the lower bound in 
the diameter of the basin theorem. □ 
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9.2 Distance to the sink (the lower bound) 

We choose a big R ^> 1 such that n = R 2 is an integer. This time we start with 

z n 

the function F(z) = —=e zRS 1 ~ rS with < 5 < 1 (later, we'll choose 5 = |). The 

Vn! 

gradient of the corresponding potential U equals 



VU (z) = ^(z)-z = --z + R 5 - 1 = ^—^ e w + R 5 - 1 , z = re w . 
F z r 

Let A = {R — 1 < \z\ < R + 1}, and let V = {z E A,\ argz — || < j^} be the same 
sector as above. Note the following properties of the gradient field: 

(i) on the boundary circumferences \z\ — R ± 1, the radial component of the field 
— Vf/ is directed outward A, and its size is not less than 1; inside A, the size 
of the radial component does not exceed 3; 

(ii) the field — Vf/ has the horizontal drift R 6 ' 1 oriented to the left; in particular, 
inside the sector V, the angular component of the field — VU is oriented counter- 
clockwise and its size is within the range 2R 5 ~ 1 ]. 

By Q we denote the set of points that hit the segment 

J= [(R- l)e l (^ + ^),(i?+l)e<^ + ^)] 

when moving along their trajectories. Because of the "left-oriented horizontal drift" 
of the field — Vf/, the points z with n ^ | arg z\ > | + ^ cannot appear within Q. (In 
fact, it is easy to see that Q C {z: < arg z < f + jq} but we will not need this). By 
G(Q) we denote the subset of Q that is located clock- wise with respect to the segment 

J (9) = [(R - l) e l ^ + To' e ) , (R + l) e <^ + TC>- e )] . 

Note that Q(9 2 ) C Q{6 1 ) for 9 2 > 9 X . By A(0) we denote the area of G(9). We denote 
by h(9) the length of the intersection of the domain Q with the segment J(0); i.e., 
the length of the "left boundary wall" of the domain Q{9). 

Lemma 9.3 If R~>\, then h{9) ^ e~ CR2 ~ S for ^ 9 ^ \, and 
m 2 (gn{-- — <&rgz <-}) ^e~ CR . 

Proof: Note that the second estimate follows from the first one by integration over 9. 
We have 

A(0) = — // AUdm 2 = — / — |<te| = — / ^-1^1 



2 .1 ]§($) 2 7^(0) (9n 2 Jg n j(0) dn 
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Figure 7: The sets Q and Q{&) 
since the rest of the boundary of Q{9) consists of gradient curves). In view of (ii) 

l -R 5 - x h{6) ^ - f 5- \ dz \ < ZR^HO) , 

2 Jar\J(9) on 



2 

whence 



ignJ(6) 

^R^htf) «C A{9) < R s ~ l h{9) . (9.4) 

We notice that \A'(9)\ ^ (R + l)h(6) < 2Rh(0). Combining this with the lower 
bound in we get the differential inequality A' (9) ^ -8R 2 ~ S A(9), whence A{9) ^ 
A(0)e- 8R2 ^ e . 

To estimate A(0) from below, recall that it equals \ the flow of the field — X7U 
through the interval J. Since the length of J is 2, A(0) cannot be less than the 
minimum of the angular component of — Vt/; i.e., A(0) ^ ^R 6 ^ 1 . Thus, A{9) ^ 

Now, using the upper bound in (19.41) . we get 

h(9) ^ R 1 - S A(9) > e- 8eR2 ~ S > e~ 2R2 - S , for ^ 9^ - . 

5 

Hence the lemma. □ 

We can replace the function F by its analytic perturbation F+H with H satisfying 
| Jr| ^ R~ 2 in the annulus R — 2 ^ \z\ ^ R + 2. After this perturbation, the gradient 
field still satisfies the conditions (i) and (ii), and the previous lemma applies to the 
new gradient flow. The next lemma gives the lower bound for the probability of the 
event that a G.E.F. / is such a perturbation. 



Lemma 9.5 If < S < 1 and R 3> 1, then P{ max 

R-2^\z\^R+2 



, I : * 
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Proof: The proof we give is very similar to that of Lemma 19. 11 Actually, we estimate 
from below the probability of the smaller event that \f — F\ ^ e~ RS \F\ everywhere in 



the annulus R — 2 < \z\ ^ R 



1 8 \z\ n 

2. Note that, in this annulus, ^ —e~ 2 

2 vn! 

dS-1 



First, we replace the exponent e zR 1 by its Taylor polynomial of degree M 
[20R 5 ] in the disk \z\ ^ 2R. It is easy to check that for m ^ M and \z\ ^ 2R, the 
m-th term in the Taylor expansion of the function e zR is bigger than twice the 
m + 1-st term. Hence the absolute value of the tail that starts with the M + 1-st 
term does not exceed the absolute value of the M-th term. In particular, the relative 
error we've made discarding the tail is at most 



Hence, for \z\ ^ 2R, 



2R s (2R 
Ml 



<S\M 



<e 2& (2eR*\ M <e -iorf 
V M J 



Pm(z)\ < e 



- 10RS \F(z) 



where 



Pm(z) 



E 

m=0 



Z \ 

m! \R, 



M 

E 



y n+m 



" y/(n + m)\ 



is the Taylor polynomial of F. Note that 

a m = R 



n + m)\ \ I _ R s 
' e 



R^ m 



The second factor on the RHS is less than 1. If R 3> 1, then the first factor does not 
exceed 2: 

) (n + l)(nt2)...M / M\ ^ 2 

V n m V #v 



Thus, < a m < 2. 
Note that 



/(*) 1 




F(z) 





f{z)-P M {z) 



F(z) 



F(z) - P w (z) 



F{z) 



and that we've already estimated the second term on the right-hand side. We write 
f(z)-P M (z) 



F(z) 



< 2e 2R max 



n<k<n+M 



E ie*-^iVnW fc ~ n + E i&i Vi 

fc^n,n+l,...,ra+M 



I fe— n 
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and show that with probability at least e CR+ ^ the maximum of the brackets on 
the right-hand side does not exceed Ce~ iR . 

We start with the first sum and demand that 

\in+m -a m \ < e- 40R \ m = 0,l, ...,M. 

The probability of this event is not less than ( ce~ 80R ) > e~ CR > e~ CR ^ . 

I Z I I x I 

For \z\ < R + 2, R 3> 1, and k > n = R 2 , we have — , ^ 2-4=. Hence 

V(fc + 1)! v^! 



fc! 



2 



1 ^ 2 fc n , and the sum we are estimating does not exceed 2 M+1 e mR * < 



e -20R" 



The second sum in the brackets does not exceed 

0^fc<n k>n+M 

We estimate from below the probability that the first sum in (19.61) is less than 3e~ mS . 
The estimate for the second sum is in the same spirit (cf. proof of Lemma I9.ip and 
we omit it. We choose a large constant A ^> 1 and split the first sum in ( 19.61) into 
two: Si = E] an d S2 = E] 

0^.k<n—A\/Mn n-A\/Mn^.k<n 

As in the proof of Lemma [9.11 we apply Lemma [2. II to estimate the sum S\. For 
this, we need to estimate the sum 

E V n ( n - 1) - (n — (n — k — 1))(R - 2) k ~ n 

0^k<n-AVMn 

E V n ( n - 1) - (n - (m - 1))(R - 2)" m . 

AV Mn<m^n 

The m-th term of the sum on the right-hand side equals 



1\ / 2\ / m-l\ ( 2 N 



1 1 ... 1 1 - - 

n J \ n J \ n J \ R 

Using inequalities 1 — £ ^ e~^, ^ £ ^ 1, and (1 — ^ e 2 ^, ^ £ ^ |, we bound 
the last expression by 

m-l 

Am -1 — r 4m (m— l)m _ m 

e r I e 2 n = e « 4 ™ ^ Ce sfl 7 ; 1 ^ m ^ n . 

j=0 
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Then the sum we are estimating does not exceed 



' / s e~^dt^CR e"* 2 dt ^ CRe~ cA2RS ^ e" 10 ^ , 



if A is big enough. 

Then, according to Lemma 12.14 the probability that 



has a double exponential decay. We conclude modestly that Si ^ e 4R with proba- 
bility at least |. 

Now, we look at the sum S2. In this case, we demand that 



|6l<e~ iaR , n - AVMn ^ k < n 
The probability of this event is not less than 

I n2(Av / M^+1) 

-■-g-lOi? 5 j > e -CR s -R 1+ ? _ e -Ci? 1+ ^ 4 



2 

Then the sum Sg does not exceed 



n-ylV Mn<fc<n 



e -io^ £ v/n(n- l)...(n- (m- l))(i2-2)' 



We know from the discussion above that each term of the latter sum is bounded by 
a constant. Hence 

if R is big enough. This completes the estimate of expression (19.61) and proves the 
lemma. □ 

Now, let us fix the variables such that the function f(z) is a small perturbation 
of F(z). For this function /, we consider the corresponding "tail" Q. If z belongs to 
the set Gn{z: § -i < argz < f } (the of this set was estimated in Lemma \9. 31) . 
then the trajectory T z must traverse the whole set Q fl {z : ~ < arg z < ~ + ^} before 
it hits the radial interval J. Hence we expect that for such z's the distance from z to 
its sink a z is comparable with R. We use this idea to prove the following lemma. 
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tt/2 + 1/10 



tt/2 - 1/10 



Figure 8: The "tail" Qi. The grey area equals Ai{9) 



i+ 3 

Lemma 9.7 Suppose R^> 1. With probability at least e~ CR 3 



— <s 

7 



m 2 {^ e D(iR, R):\z-a z \> ^} ^ e"^ 2 "* . 

Proof: After the trajectories from the tail Q leave the sector T>, they are attracted 
by some of the zeroes of the function /. Let a\, ...,ajv be the zeroes of / that 
lie in the disk 2KB and attract these trajectories, and let Qi be the corresponding 
tails. We discard the event N ^ 100i? 2 since, by Theorem 2 in (TSJ part III], its 

probability is bounded by e~ CRi which is much less than e~ CR +2 . Hence we assume 
that N < 100i? 2 . 

Let Ai(6) be the area of the tail Qi H {| — ^ < argz < | + ^ — 9}, and let ^(6*) 
be the length of the radial section of Qi by the ray { argz = | + ^ — 6*}; let A (6 I ), 
/io(#) be the similar quantities that correspond to the trajectories attracted by zeroes 
of / lying outside the disk 2RB. By Lemma 19.31 



iY 



i) = rn 2 (g n {z: f - i < argz < f }) > e" 08 " , 



=0 



thereby, Ai{jA ^ e CR for some z. 



■2-« 

10/ ^ e 

If z = 0, we are done: the points from the domain corresponding to A Q are far from 
their sinks. If i ^ 0, then, as in the proof of Lemma l9~3l hi{6) ^ R 1 ~ s Ai(9) ^ e~ CR 
for ^ 9 ^ ^j. Hence, after deleting the disk D(ai, ^), we still have a set of points 
within T> C D(iR,R) of area at least e _C7f2 * that are attracted to <2j. This proves 
the lemma. □ 
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Now, we apply the same "averaging trick" that we've already used in the proof 
of the upper bound for the distance to the sink. Consider the (random) set C = 
{z: \z - a z \ ^ ^} and the event IT = {m 2 (C fl D(z, R)) ^ e~ CR2 ~ S }. The probabil- 
ity of this event was estimated in the previous lemma (for convenience, we took there 
z = iR but, due to the translation invariance, the probability of Q* does not depend 
on the choice of z). 

We aim at estimating from below the probability P{zeC}. We have 

nR 2 F{zeC}= ¥{z + w EC} dm 2 (w) 

J J KB 

= / m 2 (CnD(z,R)) d¥^ / m 2 (C n D(z, R)) dP 
Jn Jn* 

= P { n* } e- CR2 ~ S > e -CR 1+ i s -CR^ 
It remains to put 5 — | to balance the exponents. We are done. □ 



10 Diameter of the core 

Given z G C, we show that the probability of the event {m 2 (-B z \ D(a z ,R)) > £:} 
behaves as e~ cRi when R is sufficiently large. 



10.1 The upper bound 

Given z G C, we show that the probability of the event {m2(-B 2 \ D(a z ,R)) > £:} 
cannot be bigger than e~ cRi when R ^> 1. 

We take a small positive rj depending on e only and assume that U ^ r\R 2 every- 
where in the basin B z . It is not difficult to see that the probability of the opposite 
event does not exceed Ce~ Cr,R . Indeed, the event | maxC/ > vR 2 \ is contained in 

B z J 

the union of the events { diam(S z ) > i? 4 } and { max U > r]R 2 }- By the long gra- 
dient curve theorem, the probability of the first event does not exceed Ce~ cRi . By 

2 

Lemma l4.ll the probability of the second event does not exceed CR 8 e~ ce v . 

Similarly, we also assume that U ^ —i]R 2 everywhere in B z \ D(a z , R ). The 
opposite event is contained in the union of the events { diam(S 2 ) > i? 4 } and 

{there exists a curve 7 C i? 4 D with diam(7) ^ —R such that maxt/ < —rjR 2 \ , 
and by Theorem 14. 3^ the probability of the second event is bounded by e~ Cr]Ri . 
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Thus, discarding events of probability less than e CriRi , we may assume that 
max \U\ ^ rjR 2 . Then, by our length and area estimate (Proposition! 

B z \D(a z ,R) 



{R/2f 



7]R 2 



7rexp 



□ 



m 2 {B z \D(a z ,R)) ^ 7rexp 
if 77 is sufficiently small. This proves the upper bound. 

10.2 The lower bound 

We fix a positive k < n and consider the random set C = {z : m 2 [B z \ D(a z , R)J ^ ft}. 
We need to estimate from below the probability P { z G C } , which does not depend 
on the choice of z. We apply the averaging again, but this time we average over the 
disk of radius R 5 . We get 

7iR 10 P { G C } = [[ P{ w G C } dm 2 (w) = [ m 2 (C n i? 5 D) dF . 

Introduce the event Q* that the following two conditions hold: 



(i) #(2>niflD) > 



R 2 



IX — K 



(ii) there is no gradient curve connecting the circumferences {\z\ = R} and = 
R 5 }. 

The probability of the first event is not less than e~ CRi . This estimate can be 
derived using the same techniques as in [T5l Part III] and in [TT], though it was not 
explicitly proved in these papers. To get this estimate, denote by m the least integer 



that is not less than 



Air 



7T — K 



R 2 , and estimate from below the probability that 



z m 

Sm 1 r 




> 







everywhere on the circumference {\z\ = R 2 }. We skip the estimate since it repeats 
the one used in the proof of Theorem 3 in [TT] . 

Next, by the long gradient curve theorem, the probability that the second event 
does not hold is less than e~ cR . Hence P { Q* } ^ e~ CR . 

Now, assuming that Q* happens, we can easily give a lower bound for the area 
of the set C fl R 5 3. Actually, we need to find only one basin B(a) with \a\ ^ R and 
m 2 (B(a) \ 2RB) ^ n. Then, by assumption (ii), this basin lies within the disk i? 5 D. 
Thereby, m 2 (C n R 5 B) ^ vr, and we are done: 



P{0 G C} ^ R~ w F{tt*} ^ ce 



-CR 4 
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To find a basin B(a) with \a\ ^ R and m 2 (B(a) \ 2RB) ^ «, we do a simple 
counting. Consider the basins B(a) with \a\ ^ i? but m 2 (B(a) fl 2i?B) > 7r — k. Let 
TV be the number of such basins. Comparing the areas, we get 

AttR 2 = m 2 (2RB) ^ ^m 2 (B(a) n 2KB) > (tt - k)JV; 

a 

that is, iV < 47r(7r — k)~ 1 R 2 . Hence, by assumption (i), there is at least one basin 
B(a), with \a\ ^ R and m^iB^a) \ 2KB) ^ This finishes the proof. □ 

11 Modified basins 

In this section, we prove the remaining Theorem ll.7l First, we describe a deterministic 
algorithm that "improves" partitions of the plane into domains of equal areas by 
cutting off the tentacles of the basins and re-allocating them closer to the sinks. 
Then we'll prove the probabilistic estimates for the sizes of the modified basins of our 
random partition. 

11.1 Cutting off the tentacles 

Suppose we are given a partition of the plane C = into bounded open domains 
of equal area, say tt, with marked points q, the "centers" of E^. Let 

Ri = mi{R: E % C D(c u R)}. 

Clearly, R { ^ 1. 

Given e G (0, 1), we choose the least satisfying the condition 

m 2 (Ei\D(ci,ri)) ^ 

with A = 10 4 e _1 and define the "kernef Ki — Ei H D(ci,ri) and the "tentacle" 
Ti = Ei\ D(ci, r») of the domain Ei. Note that m 2 Ti < 10~ A e. It is worth mentioning 
that this definition of the tentacle differs from the one we used in Section 8.2. Later 
on, the factor R~ 3 will help us to avoid large tangles of different tentacles. 

Proposition 11.1 Given e > 0, there exist open pairwise disjoint sets E[ with the 
following properties: 

(i) m 2 E' i = tt; 

(ii) C = E[ (up to a set of measure 0); 

i 

(Hi) m 2 (Ei n E'i) ^ tt - e; 
(iv) E[ c D( Ci ,ri + VE). 
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This proposition is useful when some of the domains Ei have long tentacles; that 
is, 7*1 <C Ri- The sets Ei may be assumed only measurable. Then the resulting sets 
E[ will be measurable too. 

Proof of Proposition [77771 " Split the plane C into standard unit squares. Suppose that 
Q is one of them. First, we check that the union of the tentacles T{ can cover only a 
small portion of the square Q: 

Lemma 11.2 

m 2 (gn(UTi)) «C — e. 

Proof: If the domain E* with Ri ^ R intersects the square Q, then Ei is contained 
in the square with side length AR + 1 homothetic to Q. Hence, comparing the areas, 
we note that 

N Q (R) = # {i : Ei n Q ± 0, Ri ^ R} ^ - (AR + l) 2 . 

7T 

Thus 

E -^ = i E ^ 

3 f°° N Q (R) 3 / , °°(4i2+l) 2 1 

= -r / , dR^ — - j-^- dR<: —e. 

A J x R 4 ^ 7iA J 1 i? 4 ^ 10 

□ 



Now, let E{ be a minimal square that is a union of several standard unit squares 
and that contains the set E^. 

Lemma 11.3 

E ^T^-e. 

i: QcEi 

Proof: Comparing the areas, we see that 

#{i:Qc Ei,Ri < i?| ^ i(4i? + 3) 2 . 

The rest is the same as in the previous lemma. □ 

Let Qi be the square that contains the center q of Ei (if q lies on the grid, it 
does not matter which one of several squares containing q to choose). For each pair 
(i, Q) with Q C Ei\ Qi, we choose a "storage" Si(Q) C Qfl (UjfTj) according to the 
following rules: 
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(a) m 2 Si(Q) = m 2 Ti] 



(b) for different z's, the storages Si(Q) are mutually disjoint; 



(c) for each pair (i,Q), the area of the storage Si(Q) is distributed between the 
kernels Kj D Q proportionally to their areas; i.e., 



m 2 (Si(Q)nK j )=m 2 (S i (Q)) 



m 2 (QnK j ) 



By Lemma 111.31 the total area within Q that we need to allocate to all the storages 
does not exceed ^, while by Lemma ft 1.21 the area of Q fl (UjKj) is not less than 
1 — Hence, we can meet the requirements (a) and (b). The requirement (c) does 
not impose any additional restriction. 

Now we describe the cut-off algorithm. It consists of countably many parallel and 
independent of each other processes. During the z-th process, for each square Q C 
the piece of the tentacle Tj HQ is re-allocated to some centers q such that KiDQ ^ 0. 
At the same time, some subsets of K\ fl Si(Q) are re-allocated to some centers c m 
whose kernels K m intersect one of the squares neighbouring Q. 

















— 



































































Qi 



































Figure 9: The square E{ and two sequences of unit squares 

We split the unit squares from E\Qi into two disjoint sequences {Qi, Q 2 , ■■■ , Q mi } 
and {Q mi+ i, Q mi + 2 , ... , Q m2 } such that in each sequence any two consecutive squares 
Qi and Qi+± have a common boundary side, and the last squares Q mi , Q m2 of each 
sequence have a common boundary side with the square Qi (see Figure 9). 

Let us call Tj \ Qi the grey area. First, for each j, 1 ^ j ^ mi, we swap the set 
Ti PI <5j with a part of the storage Si(Qj). More precisely, we 

i. choose parts of the storages G i: j C Si(Qj) with m 2 Gij = m 2 (Tj fl Q?); 

ii. decompose the tentacle T fl into disjoint union of subsets T^ij, I ^ i, with 
ra 2 T MJ = m 2 (G id nK{); 
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pieces of the kernels 
Ki and K r 



Figure 10: Putting the grey area Tj fl Qj into the storage Si(Qj) 

iii. for 1 ^ j ^ mi, re-allocate the grey area from Tj fl Qj to 

iv. for each I ^ i, remove the set K\ fl Gij from Ei, and re-allocate the set 

Ti i j of equal measure to E\. 

Now, the grey area occupies some parts of the storages Si(Qj). 

At the next step, starting with the square Qi, square after square, we move the 
grey area from Si(Qj) to S^Qj+i), until the whole grey area appears in the last 
storage Si(Q mi ) of the sequence of squares we are traversing. After that, we allocate 
the grey area to the center q. More formally, 

i. for 2 ^ j ^ mi, we choose subsets G' L j C Si(Qj) \ Gij such that raiG\ j = 

3-1 

^2m 2 G^k, and set G" d = G' Lj U Gij, G" iX = G^i; 
k=i 

ii. for 1 ^ j ^ mi — 1, we decompose the sets G"j into disjoint union of subsets 
G^ij with m 2 G ij ij = m 2 (G' i j+1 H Kj) ; 

iii. within each Ei, I ^ i, we replace the set K% fl Cij+i by the set 

l<j<mi-l 

[J °f equal measure. 

l<J<mi-l 
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Figure 11: Moving the grey area from the storages Si(Qj) to the storage Si(Q mi ) 

iv. In the end, the tentacles T^nQj, 1 ^ j ^ m 1; are cut off from Ei, and the set 
G" mi with m 2 G-' mi = 2^ m2 & H Qj) is added to E { . 

Then we apply the same process to the second sequence of squares {Q mi +i, ■ ■■ , Q m2 }- 

Note that all points re-allocated during the 2-th process will appear either in Tj 
or in one of the storages Si. Hence, due to the choice of the storages, these points are 
not displaced during the other steps. We see that for different i's the processes are 
independent of each other. 

We conclude that the new sets E[ are located in the v^-neighbourhoods of the 
kernels Ki and have the same area as E^. (a/5 is the length of the diagonal of the 
rectangle comprised of two adjacent standard squares.) By construction, 

Ei \ E\ C Tj U |J (S^nK,) . 
U,Q) 
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Due to the choice of the storages and Lemmas 111.21 and 111.3} 



U,Q) Q 



TIE 

m 2 (QnK l ) <: — 

5 



Recall that m 2 Tj ^ 10 A e. Hence, for each i, m 2 {Ei \ E[) < e. This proves Proposi- 
tion nxn □ 



11.2 Probabilistic estimate 

We fix e > and apply the cut-off algorithm to the basins B(a). The sink a is the 
"center" of B(a), R(a) is the least number R such that B(a) C D(a,R). As above, 
we set A = 10 4 e _1 , and 

r(a) = min < r : m 2 (B(a) \ D(a, r)) ^ 



AR 3 (a) 



Proposition 111 .11 gives us the modified basins B'(a) C D(a,r(a) + v5) satisfying 
conditions (i)-(iii) of Theorem II . 71 Let B'(a) be the modified basin with center at 
a that contains the origin. Since diam5'(a) ^ 2(r(a) + v5), the proof of condition 
(iv) in Theorem 11.71 boils down to the estimate 

P { r (a) > R } < e - cfi4 /(iog Rf/ 2 ( x L4 ) 

for R > A. 

Claim 11.5 P{|a| ^ |i? 4 } ^ e" cR4 . 

Proof: Assume that |a| ^ f-^ 4- Since the origin lies at the distance at most \/5 from 
the basin B(a), we know that there is a gradient curve that connects the circumfer- 
ences {\z\ = V5} and {\z\ = \R 4 }- This gradient curve connects the boundaries of 
the squares Q(0, ~^R A ) an d Q(0, 7^75 -R 4 )- By the long gradient curve theorem, the 

probability of this event is less than e~ cRi . □ 

Now, we prove (jll.4p . First, we suppose that R(ot) > R A . In view of the claim, 
we also assume that \a\ < \R A . We cover the disk D(0, \R i ) by a bounded number 
of standard squares Q(w, i? 4 ), and consider the square that contains the point 
a. We know that there is a gradient curve of diameter R(a) that terminates at the 
sink a. This gradient curve must connect dQ(w, ^= R A ) with dQ(w, ^= R 4 ). By the 

long gradient curve theorem, the probability of this event is less than e~ cRi . Hence, 
¥{R(a) > i? 4 } < e~ cR \ 
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Now, we suppose that R 4 ^ R(a). Set M = ^ i^Tr- By Lemma 14.14 throwing 

away an event of probability much less than e~ cR * , we may assume that U ^ M 
everywhere in D(0,3-R 4 ), in particular, everywhere in D(a, R(a)). By Claim [TT751 
we may assume that \a\ < i-R 4 . Hence, if min U < —M, then the disk R 4 ~B 

B(a)\D(alR) 

contains a curve of diameter at least \R where U < —M. By Theorem 14.3} the 
probability of this event does not exceed e ~ cRM3/2 , Thus, discarding the event of 
probability at most e~ /y lo s R ) ; we may assume that \U\ K M in B(a)\D(a,±R). 
Then by the length and area estimate (Proposition 18. 2p . the area of the set B(a) \ 
D(a,R) cannot exceed 

-2^ _^ D -13 



ne 2M = ttR 1J < - — — < 



AR 12 ^ AR 3 (a) ' 

provided that R ^ txA. Hence, after the events described above have been thrown 
away, we get r(a) ^ R. Therefore, the probability of the event {r(a) > R} does not 
exceed the sum of probabilities of the events thrown away, and we are done. □ 



12 Discussion and questions 

12.1 Optimal transportation to the zero set of G.E.F. 

Question 12.1 Does there exists a transportation T of the Lebesgue measure \vni 
to the random zero set Zf such that the tails sup P{ \T(z) — z\ > R} decay as e~ cR4 

as R oo ? 

Recall that the estimate e~ cR ( logfl ) 1 can be achieved by modification of the proof 
in [T5l Part II]. Note that, in view of the lower bound for the "hole probability" 
F{Z f f]RB = 0} ^ ce~ CR proved in [151 Part III], one cannot get a better estimate 
than ce . 



12.2 Length of the gradient curve and the travel time 

Given z, consider the gradient curve T z that passes through the point z. Let l z be 
the length of the part of the curve T z that starts at z and terminates at a z . 

Question 12.2 Find the order of decay of the tails P {£ z > R} as R — > oo. 

An interesting characteristic of the "random landscape" of the potential U is the 
time t z needed for the point z to roll down to the sink a z along the gradient curve T z . 
By analogy with some models from astrophysics, Michael Douglas asked us about the 
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order of decay of the tails P {r z > t} as t — > oo. Since div(VC/) = —2 everywhere on 
C\2f, one can show using Liouville's theorem that this probability equals e~ 2t (cf. 
Section 8). The length I measured along the gradient curve and the travel time r are 

connected by relation — = IVC/I. Since we know the distribution of the gradient field 

dr 



(recall that VU = — z), it looks tempting to use this information to simplify 

/ \ z ) 

the proofs of our main results and to achieve a better understanding of the properties 
of the random partition. 

12.3 Statistics of the basins 

There are several interesting questions related to the statistics of our random partition 
of the plane. We say that two basins are neighbours if they have a common gradient 
curve on the boundary. By N z we denote the number of basins B neighbouring the 
basin B z . Clearly, N z equals the number of saddle points of the potential U connected 
with the sink a z by gradient curves. Heuristically, since almost surely each saddle 
point is connected with two sinks, 

^ , mean number of saddle points per unit area 

EN Z = 2 - . 

mean number of zeroes per unit area 

Douglas, Shiffman and Zelditch proved in [3] that the mean number of saddle points 
of U per unit area is (They proved this for another closely related "elliptic model" 
of Gaussian polynomials. It seems that their proof also works for G.E.F.) Hence the 
question: 

Question 12.3 Prove thatEN z = §. 

We are also interested in the behaviour of the tails of the random variable N z : 

Question 12.4 Find the order of decay of P { N z > N } as N — ► oo. 

Another characteristic of the random partition is the number M of basins that 
meet at the same local maximum. Taking into account the result from [3] , we expect 
that its average equals 8. It is also interesting to look at the decay of the tails of M. 
Probably, some lower bound can be extracted from the analysis of perturbations of 
the polynomial z n — 1 similar to the one we did in Section 9. 

12.4 The skeleton topology 

By the skeleton of the gradient flow we mean the connected planar graph with vertices 
at local maxima of U and edges corresponding to the boundary curves of the basins. 
The graph may have multiple edges and loops. Our question is 
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Question 12.5 Are there any non-trivial topological restrictions on finite parts of 
the skeleton that hold almost surely? 

There is an interesting finite counterpart of this question. Choose N independent 
points a±, ... ,ajv uniformly distributed on the Riemann sphere C and consider the 
gradient flow of the random spherical potential 



Question 12.6 Describe all possible skeletons of the gradient flow on C of the po- 
tential V that are realized with positive probability. 
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